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I.  SUMMARY 


This  is  a  User’s  Manual  for  the  computer  program  UCIN-CABLE  III.  The  pro¬ 
gram  is  designed  to  study  the  two  dimensional  dynamics  of  submerged  and  partially 
submerged  towing  cables.  The  cables  themselves  may  have  multiple  branches. 

Figure  1  shows  a  schematic  illustration  of  a  typical  cable  system.  The  system 


can,  of  course,  be  much  simpler  than  that  shown  in  Figure  1.  For  example,  it 
could  consist  of  a  single  cable  segment  as  with  a  buoy  or  balloon  mooring  sys¬ 
tem  as  depicted  in  Figure  2. 


I.  Summary 


Buoy  or 
Balloon 


Figure  2.  A  Simple  Mooring  System. 

The  manual  provides  instruction  for  using  CABLE  to  study  such  cable  sys¬ 
tems.  It  also  provides  sample  input  and  output  data.  The  computer  program  it¬ 
self  is  written  in  FORTRAN.  Hence,  the  input  requirements  are  in  the  FORTRAN 
format. 
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II.  CAPABILITIES  OF  CABLE 


CABLE  is  designed  to  perform  the  following  dynamic  analyses:  Given:  a) 
the  physical  data  of  the  cable  (weight,  diameter,  length)  and  the  physical  data 
of  the  towed  body;  b)  the  configuration  of  the  cable  branches;  c)  the  fluid 
properties;  d)  the  motion  of  the  towed  end  of  the  cable;  and  e)  the  initial 
configuration  of  the  cable  system;  then  the  program  determines  the  kinematics 
(position,  velocity,  and  acceleration)  at  the  various  points  of  the  cable  sys¬ 
tem  as  well  as  the  cable  tension. 

CABLE  provides  the  user  with  the  following  options: 

1)  Either  English  or  metric  units  may  be  used. 

2)  The  cable  system  may  be  immersed  in  two  fluid  media  -  normally  air 

and  water. 

3)  The  fluid  media  may  be  given  a  uniform  "stream”  velocity. 

4)  Fluid  forces  may  be  applied  to  the  cable  system.  These  forces  include 

a)  Normal  drag  forces 

b)  Tangential  drag  forces 

c)  Added  mass  forces  due  to  the  cable  system  acceleration 

d)  Buoyancy  forces 

These  fluid  forces  may  be  partially  or  totally  neglected. 
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II  -  Capabilities  of  CABLE 


5)  The  towed  end  of  the  cable  may  have  arbitrary  motion  relative  to  a 

ship  frame.  The  ship  frame  itself  may  have  motion  in  a  straight  line. 


6)  Gravity  forces  may  be  included  or  neglected. 


In  addition  to  the  above  options,  the  user  may  arbitrarily  select  the  nu¬ 
merical  integration  parameters.  Certain  output  options  are  also  available. 
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III.  THEORETICAL  BASIS  OF  CABLE 


The  cable  system  is  modelled  by  a  series  of  cylindrical  links.  These 
links  are  connected  in  a  chain  to  simulate  the  cable.  Figure  3  depicts  such  a 
modelling.  The  number  of  links,  and  hence,  their  length,  is  a  user  option. 


Figure  3.  A  Finite-Segment  Model  of  a  Cable. 

The  model  of  Figure  3  is  a  "f inite-segment"  model  of  the  cable.  As  such, 
it  forms  a  "general-chain”  system  or  an  "open-tree"  system.  References  [1-3]* 
provide  an  analytical  basis  for  studying  these  systems.  The  governing  differ¬ 
ential  equations  are  obtained  using  Lagrange’s  form  of  d'Alembert’s  principle 
and  Kane’s  equations  as  exposited  by  Kane  and  others  [4-8].  This  procedure  has 
distinct  advantages  over  Newton's  laws  and  Lagrange's  equations  for  analyzing 
multibody  systems.  Specifically,  Kane’s  method  provides  for  the  automatic  elim 
ination  of  nonworking  internal  constraint  forces  while  avoiding  the  tedious 
differentiation  of  scalar  energy  functions. 


*  Numbers  in  brackets  refer  to  References  at  the  end  of  the  manual. 
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III.  Theoretical  Basis  of  CABLE 


j 

In  the  formulation  of  CABLE,  the  velocities  and  accelerations  are  computed 
through  vector  derivatives  which  may  be  evaluated  through  vector  cross  products. 
This  leads  to  algorithms  which  are  readily  converted  into  computer  subroutines. 

The  fluid  forces  on  the  CABLE  system  are  computed  using  procedures  outlined 
in  References  [9,  10,  12-15].  These  include  inertia  ("added  mass"),  drag  (nor¬ 
mal  and  tangential  for  cable  links),  and  hydrostatic  (buoyancy)  forces.  They 
are  represented  by  forces  passing  through  the  mass  centers  of  the  members  of  the 
cable  system  together  with  couples. 

Finally,  the  governing  equations  are  solved,  using  RKGS  [11],  a  fourth-order 
Runge-Kutta  numerical  integrator.  CABLE  is  written,  however,  so  that  other  nu¬ 
merical  integrators  may  also  be  used. 
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IV.  DEFINITIONS  OF  TERMS 


This  part  of  the  Manual  defines  the  terminology  used  in  CABLE.  These  defi¬ 
nitions  are  useful  in  understanding  the  input  requirements.  Underlining  is  used 
to  further  identify  terms  of  interest. 

1.  Cable  Model 


Figure  4  shows  a  finite-segment  cable  model  of  a  cable  system. 


Figure  4.  A  Finite-Segment  Model  of  a  Cable  System. 

It  consists  of  a  series  of  cylindrical  cable  links  connected  in  "tree-like" 
fashion  such  that  no  closed  loops  are  formed. 
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IV.  Definitions  of  Terms 


2 .  Cable  Link  System 

A  cable  link  system  is  a  set  of  connected  cable  links  together  with  the 
towed  bodies  as  shown  in  Figure  4.  Thus,  a  cable  link  system  is  a  multibocy 
system  as  defined  in  References  [2,  3]. 

3.  Cable  Link 


A  cable  link  is  an  individual  member  of  a  cable  link  system  (aside  from 
towed  bodies) .  A  cable  link  is  a  rigid  cylindrical  rod.  It  is  connected  to 
adjacent  links  and  the  towed  bodies  by  hinge  joints. 

4.  Reference  Link  and  Reference  Points 


The  uppermost  link  of  the  cable  link  system  is  called  the  reference  link 
as  shown  in  Figure  5.  The  upper  end  of  the  reference  link  is  called  the  system 
reference  point  Q.  In  a  typical  dynamic  configuration  of  the  system,  Q  is  given 
a  prescribed  motion  relative  to  the  mean  ship  frame.  (See  Section  7  below.) 

Also  each  of  the  links  and  the  towed  bodies  of  the  cable  link  system  has 
its  own  reference  point.  For  the  reference  link,  it  is  Q  the  system  reference 
point.  For  any  other  link,  say  L^,  it  is  the  joint  location  at  the  upper 
end  of  the  link.  For  any  towed  body,  say  B^,  it  is  the  joint  location  at 

the  cable  attach  point.  See  Figure  5. 
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IV.  Definitions  of  Terms 


Q  (System  Reference  Point) 


Reference  Link 


jed  Body 


Towed  Body, 

Figure  5.  Reference  Link,  Reference  Points,  and  Towed  Bodies. 

5 .  Towed  Bodies 

CABLE  has  provisions  for  handling  many  typical  towed  bodies  at  the  ends  of 
the  branches  of  the  cable  system.  These  bodies  fall  into  two  general  categories: 
1)  spheres;  and  2)  other  more  general  types  of  forward  towed  bodies.  CABLE  re¬ 
fers  to  the  former  as  towed  spheres,  and  the  latter  as  towed  bodies.  The  towed 
spheres  are  connected  to  the  cable  links  by  a  hinge  joint  at  the  center  of  the 
sphere.  The  towed  bodies  are  attached  to  the  cable  links  by  a  hinge  joint  at 
some  reference  point  within  the  body.  Neither  this  reference  point  nor  the  cen¬ 
ter  of  buoyancy  of  the  towed  body  need  be  located  at  its  mass  center. 
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IV »  Definitions  of  Terms 


6.  Link  Connection  Array 


The  configuration  of  the  cable  link  system  and  the  arrangement  of  the 
branches  can  be  described  by  a  link  connection  array.  To  develop  this  array, 
let  the  links  of  the  system  be  numbered  as  follows:  Let  the  reference  link  be 
link  1,  called  L^.  Next,  number  the  remaining  links,  towed  bodies,  and  towed 
spheres  in  ascending  progression  away  from  L^  through  the  branches  of  the 
system.  Finally,  let  the  mean  ship  frame  be  body  0,  called  L^ .  Figure  6 
illustrates  this  numbering  procedure  for  the  cable  system  of  Figures  4  and  5. 

O 


Figure  6.  Numbering  of  the  Cable  Link  System  of  Figures  4  and  5, 
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IV.  Definitions  of  Terms 


If  the  links,  towed  spheres,  and  towed  bodies  are  numbered  in  this  manner, 
each  link,  sphere,  and  towed  body,  except  the  reference  link,  is  connected  to 
one  and  only  one  adjacent  lower  numbered  link.  (Note,  however,  that  a  link,  for 
example,  L 2  may  be  connected  to  more  than  one  adjacent  higher  numbered  link.) 


The  array  listing  the  lower  numbered  links  is  called  the  link  connection 
array .  For  the  system  shown  in  Figure  6,  this  array  is: 

Link/Sphere/ 

Towed  Body:  1  2  3  4  5  6  7  8  9  10  11  12  13  14 

Link  Connection 

Array:  012345278  9  8  11  12  13 

Note  that  the  configuration  of  the  system  can  be  constructed  once  the  link  con¬ 
nection  array  is  known.  That  is,  there  is  an  equivalence  between  the  cable  link 
system  configuration  and  the  connection  array. 

7 •  Mean  Ship  Frame,  Inertia  Frame,  and  Coordinate  Systems 

In  CABLE,  the  reference  point  Q  at  the  upper  end  of  the  reference  link  may 
be  given  an  arbitrarily  prescribed  motion  relative  to  the  mean  ship  frame,  which 
itself  may  be  given  a  prescribed  motion  relative  to  an  inertial  reference  frame. 
This  is  depicted  in  Figure  7.  The  mean  ship  frame  is  assumed  to  be  rigidly 
attached  to  a  towing  vessel,  such  as  a  ship  or  a  helicopter. 
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IV.  Definitions  of  Terms 


Inertial  Frame  Mean  ShiP  Frame 

Figure  7.  Inertial  Frame,  Mean  Ship  Frame,  and  Cable  Link  System. 

Coordinate  Systems  are  associated  with  both  the  mean  ship  frame  and  the 
inertia]  frame,  as  shown  in  Figure  7.  In  both  frames,  the  X  directions  are  for¬ 
ward,  and  the  Z  directions  are  vertically  down.  The  motion  of  the  mean  ship 
frame  relative  to  the  inertial  frame  is  restricted  to  straight  line  motion  along 
(or  Xg).  However,  Q  may  have  arbitrary  motion  in  the  X-Z  plane  relative  to 
the  mean  ship  frame. 

Finally,  CABLE  has  the  option  of  locating  the  fluid  interface  (or  water 
level)  at  any  distance  h  below  (or  above,  if  h  is  negative)  the  mean  ship  frame. 

8.  Acceleration  Profiles 


The  specification  of  the  motion  of  reference  point  Q  relative  to  the  mean 
ship  frame  and  the  motion  of  the  mean  ship  frame  relative  to  the  inertia  frame 
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IV.  Definitions  of  Terms 


(See  Figure  7)  may  be  accomplished  either  by  using  an  acceleration  profile,  de¬ 
scribed  here,  or  by  using  a  precoded  function,  described  in  the  following  sec¬ 
tion. 


An  acceleration  profile  is  simply  a  set  of  data  representing  the  coordi¬ 
nates  of  selected  points  of  an  acceleration-time  curve.  Such  coordinates  may 
be  obtained  from  the  graph  of  the  acceleration  function. 


CABLE  has  the  capability  of  accepting  as  many  as  25  data  points  from  an 
acceleration  curve.  A  piecewise  linear  approximation  is  then  made  of  the  accel¬ 
eration  function.  For  example,  consider  the  acceleration  function  and  its  ap¬ 
proximation  shown  in  Figure  8.  The  acceleration,  velocity,  and  displacement  dur- 
ing  the  i  time  interval  are  then: 


c 
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Figure  8.  Acceleration  Profile  Approximation. 
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IV.  Definitions  of  Terms 


a  =  ai +  <r^nr  ><*  -  V 


(1) 


v  =  vi  +  a.  (t  -  t±)  +  (----  )(t  -  t±)2/2 

i+1  i 


(2) 


and 


3  ""3 

d  =  d  +  v  (t  -  t  )  +  a  (t  -  t  )2/2  +  * )  (t  -  t  )3/6 

■t  JL  X  -X-  ^  L»  «  X 

l+l  l 


(3) 


where  a^,  v^,  d^,  and  are  the  acceleration,  velocity,  displacement,  and  time 
at  the  beginning  of  the  i^  interval.  Thus,  the  entire  kinematic  profile  (dis¬ 
placement,  velocity,  and  acceleration)  is  known  when  the  are  given,  and  when 
v^  and  d^,  the  initial  velocity  and  displacement  (at  time  t^),  are  given. 

9.  Precoded  Functions 


If  the  motion  of  reference  point  Q  relative  to  the  mean  ship  frame,  or  if 
the  motion  of  the  mean  ship  frame  relative  to  the  inertia  frame  is  relatively 
simple,  it  may  be  convenient  to  describe  the  motion  using  a  precoded  function. 
For  example,  these  functions  may  be  used  to  describe  constant  or  decaying  sinu¬ 
soidal  speed  of  the  mean  ship  frame.  Specifically,  CABLE  provides  the  user  with 
the  following  precoded  function: 

f  =  f q  +  Ae^^  cos(pt  +  $)  (4) 
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where  f^,  A,  b,  p,  and  <j)  are  user  supplied  constants.  Hence,  for  a  constant 
ship  speed,  the  user  would  simply  provide  a  value  for  f  with  the  other  values 
being  zero  (which  are  the  default  values).  Specific  input  format  is  discussed 
in  the  following  part  of  the  Manual. 

10.  Fluid  Velocities 


CABLE  also  has  the  provision  of  allowing  each  of  the  fluids  to  have  a  uni¬ 
form  (constant  in  time  and  space)  ’'stream*1  velocity.  Specific  input  format  is 
described  in  the  following  part  of  the  Manual. 

11 .  Link  Orientation  Angles 

Since  the  motion  of  the  cable  system  takes  place  in  the  X-Z  plane,  the 
relative  orientation  of  two  adjacent  links  or  of  a  towed  body  and  its  adjacent 
link  can  be  described  by  a  single  angle  of  rotation  about  the  Y  axis.  Link 
orientation  angles  will  be  described  here  and  towed  body  orientation  angles  in 
the  next  section. 

Consider  a  typical  cable  link  L,  as  shown  in  Figure  9.  Let  X,  ,  Y  ,  and  Z, 

K  K  K  K 

represent  a  coordinate  frame  attached  to  L,  such  that  the  axis  of  L,  is  along 

K  K 

the  Xk  axis  as  shown.  Also  let  Y^  be  perpendicular  to  the  plane  of  motion  com¬ 
pleting  the  X^,  Y^,  Zk  right-handed  system.  (The  circular  arrow  indicates  that 
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Figure  9.  Typical  Link  and  Coordinate  Axes. 

Y k  points  out  of  the  page  -  the  direction  of  advance  of  a  right-handed  screw 

when  rotated  in  the  direction  of  the  arrow.)  Next,  consider  link  L,  and  its 

k 

adjacent  lower  link  L.  as  shown  in  Figure  10.  Let  the  axes  system  X.,  Y.,  and 

J  J  J 


H 


J 


'r 

z. 


l  r 

j 


*xi 


Figure  10.  Two  Typical  Adjoining  Links  and  Aligned  Axes. 


Z.  of  L.  be  initially  aligned  with  the  axes  of  L,  as  shown.  Now,  L,  can  be  lo- 
J  J  K  k 

cated  relative  to  L.  by  a  single  rotation  about  the  Y,  (or  Y.)  axis.  The  angle 

J  k  j 

of  rotation  is  taken  as  positive  when  the  rotation  is  right-handed  (or  dextral) 
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relative  to  the  Y  axes.  That  is,  it  is  positive  when  the  rotation  is  in  the 
direction  of  the  circular  arrows  shown.  A  positive  rotation  angle  0  is  shown 

K 

in  Figure  11 . 


Figure  11.  Positive  Link  Rotation  Angle. 

12.  Towed  Body  Orientation  Angles 

The  orientation  of  a  towed  body  relative  to  its  adjacent  lower  link  is  de¬ 
scribed  similar  to  that  for  the  cable  links.  Consider  a  typical  towed  body  as 
shown  in  Figure  12.  Let  the  X^,  Y^>  and  axes  represent  a  coordinate  frame 
located  at  the  cable  attach  point  with  the  axis  parallel  to  the  longitudinal 
axis  of  the  towed  body.  As  shown,  the  X^,  Y^,  and  axes  are  along  the  axial, 
laterial ,  and  normal  directions  of  the  towed  body. [13]  Next,  consider  towed 
body  and  its  adjacent  lower  link  L^  as  shown  in  Figure  13.  Let  the  axes 
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Figure  12.  Typical  Towed  Body  with  Coordinate  Axes. 


Figure  13.  Typical  Towed  Body  and  Adjoining  Link  with 
Aligned  Coordinate  Axes. 


system  X  ^  ,  and  of  L_.  be  initially  aligned  with  the  axes  of  B^  as  shown. 

Now,  B^  can  be  oriented  relative  to  L^.  by  a  single  rotation  about  the  (or  Y^) 
axis.  Note  that  for  typical  configurations  of  the  towed  body  ,  its  rotation 
angle  will  be  negative  as  shown  in  Figure  14. 
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Figure  14.  Typical  Towed  Body  Configuration  with 
Negative  Rotation  Angle. 

13 .  Initial  Cable  and  Towed  Body  Configuration 

The  initial  configuration  of  the  cable  system  before  the  equations  of 
motion  are  integrated  depends  upon  the  initial  values  of  the  relative  orienta¬ 
tion  angles  of  the  links  and  towed  bodies.  For  example,  if  the  initial  values 
of  all  the  orientation  angles  are  zero,  the  initial  configuration  of  a  single 
branch  cable  system  with  a  towed  body  would  be  as  shown  in  Figure  15. 
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Figure  15.  "Zero"  Initial  Configuration  of  a  Single 
Branch  Cable  with  a  Towed  Body. 


14.  Weight  Density,  Link  Length  and  Link  Diameter 


The  physical  parameters  needed  by  CABLE  to  describe  the  cable  links  are 
the  weight  density,  link  length,  and  link  diameter.  The  weight  density  is  de¬ 
fined  as  the  weight  per  unit  length. 


The  cable  links  may  have  different  lengths,  diameters,  and  weights.  That 
is,  the  cable  need  not  be  uniform  in  its  physical  properties. 


When  the  above  data  are  given,  CABLE  automatically  calculates  the  link 
mass  center  locations  and  the  link  inertia. 
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15  *  Weight  and  Diameter  of  the  Towed  Spheres 

In  addition  to  the  link  physical  parameters,  CABLE  also  requires  the  phys¬ 
ical  parameters  of  the  towed  spheres.  These  are  simply  the  sphere  weights  and 
diameters . 

16 .  Towed  Body  Physical  Data 

The  physical  data  needed  by  CABLE  to  describe  the  towed  bodies  is  somewhat 
more  extensive,  since  these  bodies  can  have  more  general  shapes.  The  towed 
body’s  weight ,  axial  length,  volume ,  inertia ,  mass  center  position  vector,  and 
buoyancy  center  position  vector  are  required.  The  body’s  inertia  is  its  mass 
moment  of  inertia  about  an  axis  perpendicular  to  the  plane  of  motion,  and  pass¬ 
ing  through  its  mass  center  CG.  The  mass  center  position  vector  is  speci¬ 
fied  by  components  of  the  position  vector  of  the  mass  center  of  the  body  rela¬ 
tive  to  the  cable  attach  point  (reference  point) ,  and  the  buoyancy  center  posi- 
tion  vector  r^  is  specified  by  components  of  the  position  vector  of  the  center 
of  buoyancy  relative  to  the  body’s  mass  center.  (See  Figure  16.)  In  both 
cases,  components  are  taken  along  the  body-fixed  coordinate  axes  with  negative 
numbers  indicating  negative  coordinate  directions. 
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Figure  16.  Typical  Towed  Body  with  Mass-Center  and 
Buoyancy-Center  Position  Vectors  and  Body-Fixed 
Coordinate  Axes. 


17.  Fluid  Forces  on  Cable  Links 


CABLE  has  the  options  of  including  up  to  four  kinds  of  fluid  forces  exerted 
on  the  cable  links:  1)  drag  forces  normal  to  the  link;  2)  drag  forces  parallel 
(tangent)  to  the  link;  3)  "added-mass"  forces;  and  4)  buoyancy  forces.  Refer¬ 
ence  [9]  describes  these  forces  in  detail. 

The  drag  forces  and  the  added  mass  forces  may  be  defined  as  follows:  Con¬ 
sider  a  typical  cable  link  as  shown  in  Figure  17.  Let  be  the  center  of  the 
link,  and  let  P  be  a  point  a  distance  x  from  as  shown.  The  added  mass  and 
drag  forces  at  P  (per  unit  length)  may  then  be  represented  as: 
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Figure  17.  Fluid  Force  at  a  Typical  Point  of  a  Cable  Link. 


£  =  a£n  +  bHxn  +  cUtUt 


(5) 


where  is  the  normal  component  of  the  fluid  velocity  relative  to  the  cable 
link  at  P,  and  the  coefficients  A,  B,  and  C  are: 


A  =  CMp(7t/4)d2 

(6) 

B  =  CNp(d/2) 

(7) 

C  =  Cxp(d/2) 

(8) 

where  p  is  the  fluid  mass  density,  and  d  is  the  diameter  of  the  cable  link.  C^, 
C^,  and  C^  are  coefficients  dependent  upon  the  Reynolds  number  of  the  fluid  flow 
past  the  cable  link.  These  coefficients  are  usually  determined  experimentally. 
In  CABLE,  the  following  values  are  used  [9,  12]: 


C 


M 


C 


N 


1.0 


0.0  for  Re  <0.1 
N  — 

0.45  +  5 . 93/ (Re.  )° ' 33  for  0.1  <  Re 
N 

1.27  for  400.  <  Re„  <  10s 

JN  ““ 

0.3  for  Re„  >  105 


N  - 


400.0 


(9) 


(10) 
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C, 


T 


0 . 0  for  Re^,  <0.1 


1.88/(ReT)°'74  for  0.1  <  ReT  <  100.55 
0.062  for  Rex  >  100.55 


(ID 


where  the  Reynolds  numbers  Re^  and  Re^  are  defined  as: 


ReN  =  pdW/M 


(12) 


and 


ReT  =  pd|vT|/M 


(13) 


where  (J  is  the  fluid  viscosity. 

In  Equation  (5),  the  first  term  is  called  the  "added-mass"  force,  the  sec¬ 
ond  term  is  called  the  "normal  drag"  force,  and  the  third  term  is  called  the 
"tangential  drag"  force.  In  CABLE,  the  system  of  added  mass  and  drag  forces 
acting  at  all  the  points  of  the  cable  link  are  replaced  by  a  single  force  pass¬ 
ing  through  together  with  a  couple.  For  additional  information,  see  Refer¬ 
ences  [9 ,  12]  . 

The  buoyancy  forces  are  due  to  hydrostatic  pressure  forces  exerted  on  the 
cable  link.  These  forces  are  normal  to  the  surfaces  of  the  link,  which  is  ex¬ 
posed  to  the  fluid.  For  typical  links,  the  resultant  hydrostatic  force  on  the 
link  is  normal  to  the  link  axis,  since  the  ends  of  the  link  are  not  exposed  to 
the  fluid.  For  these  links,  the  system  of  hydrostatic  forces  may  be  replaced 
by  a  single  force  normal  to  the  link  axic  given  by: 


Bn  =  PgV^  *  Qc  x  *0 


(14) 


24 


IV.  Definitions  of  Terms 


where  p  is  the  fluid  mass  density,  g  is  the  gravity  constant,  V,  is  the  cable 

K 

link  volume,  n  is  a  unit  vector  along  the  link  axis,  and  £  is  a  vertical  unit 
vector  as  shown  in  Figure  17.  CABLE  treats  all  cable  links  as  suggested  above 
except  those  at  the  end  of  a  branch  which  has  no  towed  sphere  or  towed  body. 

For  those  segments  an  adjustment  is  made  to  the  buoyancy  force  to  include  the 
hydrostatic  forces  on  the  exposed  end. 

(For  additional  information,  see  reference  [9]) 

18 .  Fluid  Forces  on  Towed  Spheres 

CABLE  has  the  option  of  including  up  to  three  kinds  of  fluid  forces  exerted- 
on  the  spheres:  1)  drag  forces;  2)  "added-mass"  forces;  and  3)  buoyancy  forces. 
CABLE  replaces  the  system  of  drag  forces,  added  mass  forces,  and  buoyancy  forces 
by  a  single  force  passing  through  the  mass  center  of  the  sphere. 

The  resultant  drag  force  and  added  mass  force  is  taken  to  be  [12,  15]: 

F  =  -Aa  +  b|v|v  (15) 

where  a  is  the  acceleration  of  the  mass  center  of  the  sphere  relative  to  the 
fluid,  V  is  the  fluid  velocity  relative  to  the  mass  center  of  the  sphere,  and 
the  coefficients  A  and  B  are: 
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A  = 

(16) 

B  =  Vd2CD 

(17) 

Here  p  is  the  fluid  mass  density,  is  the  volume  of  the  sphere,  and  d  is  the 
diameter  of  the  sphere.  and  are  coefficients  dependent  on  the  Reynold's 

number  of  the  fluid  flow  past  the  sphere.  These  coefficients  are  usually  deter¬ 
mined  experimentally.  In  CABLE  the  following  values  are  used  [12]: 


C 


M 


1.0 


0.0  for  Re  <  0.1 


CD  = 


0.044  +  13 . 46/ (Re) 2  for  0.1  <  Re  <  1000.0 

0.47  for  1000.0  <  Re  <  105 
0.12  for  Re  <  105 


where  the  Reynold's  number  Re  is  defined  as 
Re  =  pd|v|/M 

where  p  is  the  fluid  viscosity. 


(18) 


(19) 


(20) 
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The  buoyancy  forces  are  due  to  the  hydrostatic  pressure  forces  exerted  on 
the  sphere.  These  forces  are  normal  to  the  surface  of  the  sphere  which  is  ex¬ 
posed  to  the  fluid.  Since  CABLE  assumes  the  cable  links  to  be  imbedded  directly 
into  the  sphere,  it  modifies  the  usual  buoyancy  force  on  the  sphere  to  account 
for  the  small  portion  of  its  surface  which  is  shielded  from  the  fluid.  In 
general,  then,  the  resultant  buoyancy  force  is  not  directed  exactly  vertical. 

19.  Fluid  Forces  on  the  Towed  Bodies 


CABLE  has  the  option  of  including  up  to  three  kinds  of  fluid  forces  on  the 
towed  bodies:  1)  drag  forces;  2)  "added-mass”  forces;  and  3)  buoyancy  forces. 
The  drag  and  added  mass  forces  are  calculated  from  a  set  of  25  hydrodynamic  co¬ 
efficients.  The  forms  of  these  forces  were  developed  from  References  [13,  14]. 

Before  defining  the  drag  and  added  mass  forces,  first  consider  the  typical 
<owed  body  shown  in  Figure  18. 
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0 


s 


X 

(forward) 


u»u 


Z 


Figure  18.  Typical  Towed  Body  with  Positive  Coordinate, 
Velocity,  Angular  Velocity,  and  Angle  Directions. 


Here  u  and  w  represent  the  velocity  components  of  the  mass  center  CG  of  the 
body  along  the  body-fixed  directions,  q  represents  the  angular  velocity  of  the 


body,  and  6^  and  represent  the  angle  deflections  of  the  bowplane  and  stern- 


plane,  respectively.  The  directions  associated  with  the  body-fixed  X,  Y,  and 
Z  coordinate  axes  are  called  the  axial ,  lateral ,  and  normal  directions. 

CABLE  replaces  the  drag  and  added  mass  forces  by  a  single  force  and  a 
couple  which  are  defined  as  follows: 

Axial  Force 


(21) 


+  hp  &2  [X'  u2  +  X'  w2  +  Xl  u262  +  X~  u262] 
r  uu  ww  6,  B  o  s 


b  s 
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Normal  Force 


Fz  -  hp  a4  [Z^q]  +  hp  ^3[Z^w  +  Zjquq  +  Z  'q  Jvj^q  |w|  ] 

+  hp  £2  [Z'  u2  +  Z'  uw  +  Z'  u26  +  Z'  u26  +  Z"  fw|w|  ] 
uu  uw  6^  b  6g  s  w|w|  1  1  J 


(22) 


Pitch  (Lateral)  Moment 

M^  =  hp  &5  [M*q]  +  hp  &4  [Mrw  +  M"  uq  +  M',  ,q|w|  ] 
i  q  w  uq  ^  q  |  w|  ^ ' 

+  hp  &3  [M"  u2  +  M'  uw  +  M^  u26,  +  M"  u26 
uu  uw  6,  b  6  s 

b  s 

+  M  '  ,  w  |w|  ] 

w  |w| 


(23) 


where  p  is  the  mass  density  of  the  fluid  and  £  is  the  axial  length  of  the  towed 
body.  The  X  ,  Z  ,  and  M'  with  the  appropriate  subscripts  are  the  hydrodynamic 
coefficients  for  the  axial,  normal,  and  lateral  directions,  respectively. 

(These  coefficients  must  be  supplied  by  the  user  on  input.)  It  should  be  noted 
here  that  the  body  has  been  assumed  symmetric  with  respect  to  the  X-Z  plane  and 
that  it  is  undergoing  a  fully  submerged  forward  tow. 


The  bowplane  and  sternplane  deflections ,  6^  and  6g ,  must  also  be  specified 
by  the  user.  These  are  not,  however,  "read  in"  with  the  normal  input  data,  but 
are  specified  in  SUBROUTINE  CVSPEC.  The  subroutine  supplied  with  CABLE  simply 
sets  the  values  of  6^  and  6 to  zero  and  returns  to  the  calling  segment.  This 
subroutine  may  be  easily  replaced  to  give  the  deflections  nonzero  values. 
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The  subroutine  statement  must  read: 

SUBROUTINE  CVSPEC  (TIME,  I BODY,  IFL,  U,  W,  Q,  THETA,  X,  Z,  DELS,  DELB) 

where  the  variables  in  the  parameter  list  are  defined  as  follows:  (*  -  indi¬ 
cates  the  value  is  supplied  by  CABLE) 


VARIABLE 

MEANING 

TIME* 

Value  of  time  in  seconds 

I BODY* 

Integer  body  number  of  the  towed  body 

IFL* 

Integer  specifying  the  towed  body  to  be  in 
Fluid  1  (IFL  =  1)  or  Fluid  2  (IFL  =  2). 

u* 

Axial  component  of  the  velocity  of  the  mass- 
center  of  the  towed  body. 

W* 

Downward  (Z)  component  of  the  velocity  of 
the  mass-center  of  the  towed  body 

Q* 

Angular  velocity  of  the  towed  body 

THETA* 

The  angle  between  the  towed  body  axes  and 
the  inertial  axes  (positive  rotation  in  the 
+Y) . 

X* 

The  X  coordinate  of  the  mass-center  relative 
to  the  inertial  frame. 

z* 

The  Z  coordinate  of  the  mass-center  relative 
to  the  inertial  frame. 

DELS 

Sternplane  deflection  (radians)  (See  Fig¬ 
ure  18.) 

DELB 

Bowplane  deflection  (radians)  (See  Fig¬ 
ure  18.) 
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The  values  of  DELS  (6  )  and  DELB  (6_)  must  be  calculated  and  returned  to  CABLE. 

S  D 

Finally,  since  CABLE  assumes  that  the  center  of  buoyancy  and  the  center  of 
mass  of  the  towed  body  do  not  necessarily  coincide,  it  replaces  the  buoyancy 
forces  by  a  force  through  the  mass  center  of  the  body  as  well  as  a  couple. 

20.  Weight  Forces 


The  weight  force  may  be  represented  on  each  cable  link  L,  by  a  single 
vertical  (downward)  force  passing  through  the  mass  center  of  If 

is  the  weight  density  per  unit  length  of  Lt  ,  W,  may  be  expressed  as: 

K 

4  =  -W  <24> 

where  i-s  the  length  of  and  where  ^  is  the  vertical  unit  vector  shown  in 
Figure  17. 


Similarly,  the  weight  forces  on  the  towed  spheres  and  towed  bodies  may  be 
represented  by  single  vertical  forces  passing  through  the  center  of  the  bodies 
with  magnitude  equal  to  their  weights. 

21.  Units 


CABLE  allows  the  user  to  select  either  English  or  metric  units  for  the  in¬ 
put  and  output  data.  If  English  units  are  selected,  the  units  are  slugs ,  feet , 


31 


IV.  Definition  of  Terms 


and  seconds  for  mass,  length,  and  time.  However,  inches  are  used  for  the  link 
diameter  and  towed  sphere  diameter.  If  metric  units  are  selected  the  units  are 
kilograms ,  meters ,  and  seconds  with  centimeters  used  for  link  diameter  and  towed 
sphere  diameter. 

22.  Labels 


CABLE  allows  the  user  to  arbitrarily  label  or  name  the  links  and  towed 
bodies,  their  reference  points,  and  the  towed  spheres. 
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This  part  of  the  Manual  describes  the  specific  input  data  requirements  of 
CABLE.  The  data  itself  is  described  in  terms  of  "line  images"  such  as  would  al¬ 
so  appear  on  a  computer  card. 

As  noted  earlier,  the  algorithms  and  code  of  CABLE  are  written  in  FORTRAN. 
Hence,  the  input  data  is  to  be  submitted  in  the  FORTRAN  format.  Integer  data 
is  to  be  entered  in  15  format,  and  real  data  in  either  F10.3  or  E10.3,  unless 
otherwise  specified.  For  example,  if  two  real  numbers  are  to  be  read  on  a 
given  line,  then  2F10.3  or  2E10.3  will  be  used.  The  order  of  the  input  data 
is  in  the  same  order  as  it  is  described  below.  The  units  of  the  input  data  are 
either  English  or  metric,  according  to  the  option  selected. 

1 .  Heading  or  Title 


The  first  line  of  input  to  CABLE  contains  a  title  which  may  be  used  to 
identify  the  particular  data  entered.  It  may  contain  up  to  80  characters,  and 
it  will  appear  at  the  top  of  each  computer  output  data  page. 

2 .  Options 

CABLE  has  a  number  of  user  options  to  define  the  desired  analysis.  These 
options  are  selected  by  specifying  integer  values  (usually  0,  1,  2,  3,  4,  or  5) 
as  follows: 
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V.  Input  Data 


1)  Units  (See  Section  21.,  Part  IV):  On  the  first  line  of  the 
data,  enter  a  1  if  English  units  are  desired.  Enter  a  2  if 
metric  units  are  desired. 

2)  Normal  drag  forces  (See  Sections  17,  18,  19.,  Part  IV):  Enter 
a  1  to  include  these  forces.  Enter  a  0  to  neglect  these  forces 

3)  Tangential  drag  forces  (See  Sections  17,  18,  19.,  Part  IV): 
Enter  a  1  to  include  these  forces.  Enter  a  0  to  neglect  these 
forces . 

4)  Added  mass  forces  (See  Sections  17,  18,  19.,  Part  IV):  Enter 
a  1  to  include  these  forces.  Enter  a  0  to  neglect  these  forces 

5)  Bouyancy  forces  (See  Sections  17,  18,  19.,  Part  IV):  Enter  a 
1  to  include  these  forces.  Enter  a  0  to  neglect  these  forces. 

6)  Gravity  forces  (See  Section  20.,  Part  IV):  Enter  a  1  to  in¬ 
clude  these  forces.  Enter  a  0  to  neglect  these  forces. 

7)  Viscosities  and  fluid  mass  densities  (See  Sections  17,  18,  19., 
Part  IV):  Enter  a  0  to  use  the  default  values  coded  in  CABLE. 
See  Table  I.  If  values  different  from  these  default 
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V.  Input  Data 


values  are  desired,  enter  a  1.  Then  four  lines  of  data  must 
be  entered,  providing  the  following: 

i)  Viscosity  of  the  first  fluid  (usually  air)  (lb. sec. /ft2 
or  N. sec . /m2) . 

ii)  Mass  density  of  the  first  fluid  (slug/ft3  or  kg/m3). 


TABLE  I.  Default  Values  for  Fluid  Viscosities 
and  Mass  Densities 

English  Metric 


Viscosity  (air) 

3.7188  x  10“7(lb-sec)/ft2 

1.7802  x  10~s(N-sec)/m2 

Mass  (air) 

Density 

2.378  x  10~3  slug/ ft3 

1.2252  kg/m3 

Viscosity  (water) 

3.516  x  10~5(lb-sec)/ft2 

1.6831  x  10  3(N-sec)/m2 

Mass  (water) 

Density 

1.9856  slug/ft3 

1.0231  x  103  kg/m3 

iii)  Viscosity 

of  the  second  fluid  (usually  water)  (lb. sec. /ft2 

or  N. sec/m2) . 


iv)  Mass  density  of  the  second  fluid  (slug/ft3  or  kg/m3). 


8)  Fluid  velocities:  CABLE  has  the  option  of  allowing  the  fluids 
to  have  constant  stream  velocities  in  any  direction.  To  use 


35 


V.  Input  Data 


this  option,  enter  a  1.  Next,  enter  two  lines  of  data  provid¬ 
ing:  the  X  and  Z  components  relative  to  inertia  space  of  the 

velocities  of  the  two  fluids  (fluid  1  (air)  on  the  first  line, 
and  fluid  2  (water)  on  the  second  line) . 

To  decline  this  option,  enter  a  0.  This  means  the  fluids  will 
have  zero  velocity. 

9)  Fluid  Interface  Level:  CABLE  has  the  option  of  allowing  the 
fluid  interface  level  (that  is,  the  water  level)  to  be  located 
at  some  nonzero  vertical  distance  below  the  mean  ship  frame. 

To  use  this  option,  enter  a  1.  Then,  on  the  next  line  of  in¬ 
put,  enter  the  vertical  distance  (positive  downward)  from  the 
mean  ship  frame.  To  decline  this  option,  enter  a  0.  This 
means  the  fluid  interface  and  the  mean  ship  frame  will  be  on 
the  same  level. 

The  set  of  options  data  will  consume  between  9  and  16  lines ,  depending  upon 
whether  viscosities,  densities,  velocities,  and  the  fluid  interface  level  aiS 
specified  by  the  user* 
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V.  Input  Data 


3.  Number  of  Links 


Following  the  input  on  user  options,  enter  the  number  of  cable  links  on  a 
single  line. 

4.  Individual  Link  Data 


The  next  input  data  is  a  series  of  data  sets  describing  the  physical  and 
geometrical  properties  of  each  link,  as  follows: 


1)  Link  number  (See  Section  6.,  Part  IV). 


2)  Label:  Any  name  up  to  20  characters  in  length. 


3)  Link  Reference  Point  Label  (or  "joint"  label)  (See  Section  4., 
Part  IV):  Any  name  up  to  20  characters  in  length. 


4)  Link  number  of  adjacent  lower  numbered  link. 


5)  Link  weight  density  (weight  per  unit  length,  lb/ft  or  N/m). 


6)  Link  length  (ft  or  m) . 
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V.  Input  Data 


7)  Link  diameter  (in  or  cm). 

8)  Initial  configuration  of  the  link.  This  line  contains  a  single 
number  representing  0,  the  angle  of  the  link  (in  degrees)  rela- 
tive  to  its  adjacent  lower  numbered  link  (See  Sections  11.  and 
13. ,  Part  IV). 

9)  Initial  motion  of  the  link.  This  line  contains  a  single  number 
representing  0,  the  derivative  of  the  initial  angle  of  the  link 
(in  degrees/second). 

Nine  lines  of  data  are  required  for  each  link. 

5 .  Number  of  Towed  Spheres 

The  link  data  is  followed  by  analogous  data  for  the  towed  spheres.  The 
first  line  of  towed  sphere  data  contains  the  number  of  towed  spheres. 

6 .  Individual  Towed  Sphere  Data 

The  next  input  data  is  a  series  of  data  sets  describing  the  physical  and 
geometrical  properties  of  each  towed  sphere  as  follows: 
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V-  Input  Data 


1)  Towed  sphere  number  (See  Section  6.,  Part  IV). 

2)  Towed  sphere  label  (up  to  20  characters). 

3)  Link  number  of  the  adjacent  towing  link. 

4)  Towed  sphere  weight  (lb  or  N) . 

5)  Towed  sphere  diameter  (in  or  cm). 

Five  lines  of  data  are  required  for  each  towed  sphere. 

7 •  Number  of  Towed  Bodies 

The  towed  sphere  data  is  followed  by  analogous  data  for  the  towed  bodies. 
The  first  line  of  the  towed  body  data  contains  the  number  of  towed  bodies. 

8*  Individual  Towed  Body  Data 

The  next  input  data  is  a  series  of  data  sets  describing  the  physical,  geo^ 
metric,  and  hydrodynamic  properties  of  each  towed  body  as  follows: 

1)  Towed  body  number  (See  Section  6.,  Part  IV). 
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V.  Input  Data 


2)  Label:  Any  name  up  to  20  characters  in  length. 

3)  Towed  body  reference  point  label  (or  "joint"  label)  (See  Sec¬ 
tions  4.,  5.,  and  16.,  Part  IV):  Any  name  up  to  20  characters 
in  length. 

4)  Link  number  of  adjacent  towing  link  (See  Section  6.,  Part  IV). 

5)  Towed  body  weight  (lb  or  N) . 

6)  Axial  length  of  body  (ft  or  m) . 

7)  Towed  body  volume  (ft3  or  m3). 

8)  Towed  body  inertia  (See  Section  16.,  Part  IV):  (slug-ft2  or 
N-m2) . 

9)  Mass  center  position  vector.  This  line  contains  two  numbers 
representing  the  body-fixed  X  and  2  components  of  the  position 
vector  of  the  mass  center  of  the  towed  body  relative  to  the 
cable  attach  point  (ft  or  m) .  Negative  numbers  indicate  nega¬ 
tive  coordinate  directions.  (See  Section  16.,  Part  IV). 
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V.  Input  Data 


10)  Center  of  buoyancy  position  vector.  This  line  contains  two 
numbers  representing  the  body- fixed  X  and  Z  components  of  the 
position  vector  of  the  center  of  buoyancy  of  the  towed  body 
relative  to  its  mass  center  (ft  or  m) .  Negative  numbers  indi¬ 
cate  negative  coordinate  directions.  (See  Section  16.,  Part  IV). 

11)  Initial  configuration  of  the  towed  body.  This  line  contains  a 
single  number  representing  0,  the  initial  angle  of  the  towed 
body  (in  degrees)  relative  to  its  adjacent  towing  link.  (See 
Sections  12.  and  13.,  Part  IV). 


12)  Initial  motion  of  the  towed  body.  This  line  contains  a  single 
number  representing  0,  the  derivative  of  the  initial  angle  of 
the  towed  body  (in  degrees/second). 


13)  Hydrodynamic  coefficients  for  the  towed  body.  The  next  six 
lines  of  input  contain  the  values  representing  the  hydrody¬ 
namic  coefficients  used  to  calculate  the  axial  and  normal 
forces  and  the  pitch  (lateral)  moment  on  the  towed  body. 
They  appear  on  the  lines  as  follows: 


Line  1 : 


X" 


ww 


and  X' 

qq 


Line  2: 


and  Xg 

s 
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V.  Input  Data 


Line  3:  ZT 


w 


Line  4: 


s 


Line  6: 


s 


(See  Section  19.,  Part  IV). 


Eighteen  lines  of  data  are  required  for  each  towed  body. 

9 .  Motion  of  the  Mean  Ship  Frame  and  the  System  Reference  Point 

The  mean  ship  frame  may  have  straight  line  motion  relative  to  the  inertial 
frame  R.  One  variable  representing  the  forward  (or  reverse)  speed  are  used  to 
describe  the  motion.  Also,  the  system  reference  point  Q  may  have  arbitrary  mo¬ 
tion  relative  to  the  mean  ship  frame,  requiring  two  additional  variables.  (See 
Section  7.,  Part  IV). 


For  input,  these  variables  are  identified  by  the  following  integers: 
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V.  Input  Data 


1  -  Forward  displacement  (X)  of  the  reference  point  Q  relative  to 

the  mean  ship  frame  (backward  if  negative). 

2  -  Vertical  (downward)  displacement  (Z)  of  the  reference  point 

point  Q  relative  to  the  mean  ship  frame. 

3  -  Forward  (or  reverse,  if  negative)  speed  of  the  mean  ship  frame  rela¬ 

tive  to  the  inertial  frame  R.  Note:  If  a  precoded  function  is  used 
to  specify  this  motion,  the  mean  ship  frame  is  assumed  to  coincide 
with  the  inertial  frame  at  t  =  0.0. 

The  data  describing  these  variables  may  be  specified  in  two  ways:  1)  by 
using  precoded  functions  and  2)  by  using  acceleration  profiles.  (See  Sections 
8.  and  9.,  Part  IV.)  The  precoded  functions  are  of  the  form: 

f  =  fQ  +  Aebtcos(pt  +  <j))  (25) 

where  f^  +  A,  b,  p,  and  <))  are  constants.  The  acceleration  profiles  are 
described  in  Section  8.  of  Part  IV  of  the  Manual. 

The  specific  input  data  requirements  are  as  follows:  First,  enter  the 
number  of  these  variables  which  are  to  have  nonzero  values  during  the  motion. 
Then,  for  each  of  these  variables  enter  its  identifying  integer  (see  above) 
in  column  1.  In  column  2  enter  a  1  or  a  2  depending  on  whether  a  precoded 
function  or  an  acceleration  profile  is  desired. 
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V.  Input  Data 


If  a  precoded  function  is  selected,  enter  the  values  of  f^,  A,  b,  p,  and 
<{>  on  the  next  line. 

Alternatively,  if  an  acceleration  profile  is  selected,  enter  the  number 
of  data  points  (up  to  25)  on  the  next  line.  This  line,  in  turn,  must  then  be 
followed  by  the  same  number  of  data  lines  as  there  are  data  points  (one  line 
per  point).  The  first  of  these  lines  contains:  i)  the  initial  time,  ii)  the 
initial  acceleration,  iii)  the  initial  speed,  and  iv)  the  initial  displace¬ 
ment.  The  subsequent  lines  contain  two  numbers:  i)  the  time  and  ii)  the 
acceleration  value. 

If  several  variables  have  specified  motion,  their  data  is  included  in 
serial  fashion. 

10 .  Integration  Options 

The  next  data  is  a  set  of  four  numbers  setting  the  parameters  of  the 
integration  subroutine  RKGS.  These  numbers  are  entered  on  a  single  line  as 
follows:  i)  Integration  starting  time,  ii)  Integration  ending  time,  iii) 

Integration  time  increment,  and  iv)  Upper  bound  on  the  error.  Note:  If  the 
integration  time  increment  is  input  as  a  negative  value,  say  -I,  then  CABLE 
will  automatically  reset  the  increment  to  (2)  ^ . 
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V,  Input  Data 


11 .  Output  Options 

1)  The  next  line  contains  the  time  increment  for  printing  the  computed 
data.  Since  printing  this  data  at  each  integration  step  can  produce 
an  excessive  amount  of  printed  output,  the  print  time  increment  is 
usually  greater  than  the  integration  time  increment. 

2)  At  the  beginning  of  the  output  data,  CABLE  lists  a  copy  (or  "echo") 
of  the  input  data.  The  extent  of  the  listing  of  computed  data  is 

a  user  option.  To  exercise  this  option,  enter  an  integer  (0,1, 2, 3, 4) 
in  column  1  according  to  the  following  printing  priority: 

0  -  Prints  all  compute!  data:  Position,  velocity,  and  acceleration 
of  the  mean  ship  frame;  Position,  velocity,  and  acceleration  of  Q 
the  system  reference  point;  Orientation  angles  and  their  deriva¬ 
tives  for  each  link  and  towed  body  (in  degrees,  degrees/second ,  and 
degrees/second/second);  Position,  velocity,  and  acceleration  of  the 
link,  towed  sphere,  and  towed  body  mass  centers  relative  to  inertia 
space;  Position  velocity,  and  acceleration  of  the  connecting  joints 
relative  to  the  inertia  space;  Cable  tension  at  the  system  reference 
point  Q;  and  Cable  tension  at  the  connecting  joints. 
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V.  Input  Data 

1  -  Prints  as  above  for  the  option  integer  0,  except  for  the  cable  tension 

at  the  connection  joints. 

2  -  Prints  as  above  for  the  option  integer  1,  except  for  the  cable  tension 

at  the  system  reference  point. 

3  -  Prints  as  above  for  the  option  integer  2,  except  for  the  position 

velocity,  and  acceleration  of  the  joints. 

4  -  Prints  as  above  for  the  option  integer  3,  except  for  the  position, 

velocity,  and  acceleration  of  the  mass  centers. 
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VI.  EXAMPLE  INPUT  DATA 


For  an  illustration  of  the  input  data  consider  the  cable  link  system 
shown  in  Figure  19. 


Let  the  eleven  cable  links  be  identical  having  length  2  ft,  diameter  0.25  in., 
and  weight  density  of  0.2  lb/ft.  Let  the  two  towed  spheres  (Si  and  S2)  have 
diameters  1.0  in.  and  1.5  in.  with  weights  of  0.5  lb  and  1.688  lb.  Let  the 
towed  body  (B^  be  a  sphere  of  diameter  2.0  in.  (0.166  ft),  volume  0.0024241 
ft3,  weight  4  lb  and  inertia  3.4507  x  10  4  slug-ft2.  Finally  let  the  reference 
point  be  on  the  surface  of  the  sphere,  and  let  the  centers  of  mass  and  buoyancy 
both  be  located  at  the  center  of  the  sphere. 


47 


VI.  Example  Input  Data 


Let  the  mean  ship  frame  have  a  forward  motion  defined  by  the  acceleration 
profile  shown  in  Figure  20.  Let  the  initial  forward  ship  speed  be  5  ft/sec, 
and  let  the  initial  displacement  be  zero. 


Time  (sec) 

Figure  20.  Acceleration  Profile  for  Forward  Mean  Ship  Frame  Motion 

Let  the  system  reference  point  (point  A)  have  oscillatory  motion  relative 
to  the  mean  ship  frame  given  by  the  expression 

z  =  0.5  cos  Ttt 

Let  the  fluid  level  be  3  ft  below  the  mean  ship  frame. 

Let  the  links,  spheres,  and  towed  bodies  be  labelled  as  shown  in  Figure  19 
and  numbered  as  shown  in  Figure  21. 
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Let  the  system  be  moving  in  still  air  and  calm  water.  Let  gravity  and 
fluid  forces  be  acting  and  let  the  air  and  water  viscosities  and  densities  be 
the  same  as  the  default  values  in  CABLE. 

Finally,  to  illustrate  the  input  of  initial  conditions,  let  link  make 
an  angle  of  45°  with  the  horizontal  and  let  link  L2  make  an  angle  of  30°  with 
Lx  as  shown  in  Figure  19.  Let  the  initial  angles  of  all  the  other  links  be 
zero  (even  though  they  are  not  depicted  that  way  in  Figure  19.) 

The  input  data  for  the  two-dimensional  analysis  of  the  motion  is  then  as 
follows:  (Note:  The  comments  in  parentheses  are  for  illustration  purposes 

only  and  should  not  appear  in  the  normal  FORTRAN  input.) 
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VI.  Example  Input  Data 


Line  Column  Numbers 


12  3 

123456789012345678901234567890123. .  . 

EXAMPLE  INPUT  DATA 
1 
1 
1 
1 
1 
1 
0 
0 
1 

3.0 

11 

1 

LI 

A 

0 

0.2 

2.0 

0.25 

45.0 

0.0 

2 

L2 

B 

1 

0.2 

2.0 

0.25 

-30.0 

0.0 

3 

L3 

C 

2 

0.2 

2.0 

0.25 

0.0 

0.0 

12  3 

123456789012345678901234. .  .  0 


(Title) 

(English  Units) 

(include  normal  drag) 

(include  tangential  drag) 
(include  added  mass  effects) 
(include  buoyancy  forces) 
(include  gravity  forces) 

(use  default  fluid  properties) 
(fluids  are  at  rest) 

(specify  by  nonzero  fluid 
interface  level) 

(number  of  cable  links) 

(link  body  number) 

(link  label) 

(reference  point  label) 

(lower  link) 

(weight/ length) 

(length) 

(diameter) 

(initial  orientation) 

(initial  angle  derivatives) 

(next  link  body  number  and  data) 


(next  link  body  number  and  data) 
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VI .  Example  Input  Data 


Line  Column  Numbers 

1  2  3 
123456789012345678901234567890123. 

4 

L4 

D 

3 

0.2 

2.0 

0.25 

0.0 

0.0 

5 

L5 

E 

4 

0.2 

2.0 

0.25 

0.0 

0.0 

7 

L6 

C 

2 

0.2 

2.0 

0.25 

0.0 

0.0 

8 

L7 

F 

7 

0.2 

2.0 

0.25 

0.0 

0.0 

9 

L8 

G 


(next  link  body  number  and  data) 


(next  link  body  number  and  data) 


(next  link  body  number  and  data) 


(next  link  body  number  and  data) 


(next  link  body  number  and  data) 


12  3 

123456789012345678901234...  0 
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VI .  Example  Input  Data 

Line  Column  Numbers 

12  3 

123456789012345678901234567890123. 
8 

0.2 

2.0 

0.25 

0.0 

0.0 

11 

L9 

G 

8 

0.2 

2.0 

0.25 

0.0 

0.0 

12 

L10 

H 

11 

0.2 

2.0 

0.25 

0.0 

0.0 

13 

Ill 

I 

12 

0.2 

2.0 

0.25 

0.0 

0.0 


2 

6 


51 

0.5 

1.0 

52 


10 


12  3 

123456789012345678901234...  0 


(next  link  body  number  and  data) 


(next  link  body  number  and  data) 


(next  link  body  number  and  data) 


(number  of  towed  spheres) 
(sphere  body  number) 

(label) 

(body  number  of  lower  link) 
(weight) 

(diameter) 

(next  towed  sphere  body  number 
and  data) 
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VI .  Example  Input  Data 


Line  Column  Numbers 

12  3  4 

1234567890123456789012345678901234567890123. . . 


1.688 

1.5 


1 


14 

B1 

J 

13 

4.0 

0 . 1666667 
0.0024241 

3.4507E- 

4 

-0.0833333  0.0 

0.0 

0.0 

0.0 

0.0 

-0.52 

-0.46 

-0.23 

0.0 

0.0 

0.0 

-0.52 

0.0 

0.0 

-0.46 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

21 

0.0 

0.5 

0.0 

3.1416 

32 

0.0 

0.0 

5.0 

0.0 

5.0 

2.0 

15.0 

0.0 

0.0 

20.0 

0.25 

0.001 

0.5 

0 

1 

2 

3 

123456789012345678901234. . . 


(number  of  towed  bodies) 
(towed-body  body  number) 

(body  label) 

(joint  label) 

(lower  link  body  number) 

(weight) 

(axial  body  length) 

(volume) 

(inertia) 

(c.m.  position  vector) 

(c.b.  position  vector) 

(initial  angle) 

J  (initial  angle  derivative) 

(hydrodynamic  coefficients) 

(number  of  specified  variables) 
(variable  2  is  specified  with  a 
0.0  precoded  function) 

(variable  3  is  specified  with  an 
acceleration  profile) 

(integration  parameters) 

(time  increment  for  printing  data) 
(output  option--print  all  data) 
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VII.  EXAMPLE  OUTPUT 


Besides  providing  an  "echo"  of  the  input  data  prior  to  numerical  integra¬ 
tion,  CABLE  provides  output  data  at  specified  instants  of  time.  Both  the  print¬ 
ing  interval  size  as  well  as  the  amount  of  printout  received  at  the  end  of  the 
interval  are  user  options.  Shown  in  Figures  22  and  23  is  the  standard  printout 
for  an  output  printing  priority  of  zero. 
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Figure  23.  Example  Output 


APPENDIX 


SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


This  Appendix  to  the  Manual  provides  an  abbreviated  summary  of  the  input 
data  requirements  for  a  quick  reference. 


1*  Heading  or  Title:  Heading  of  up  to  80  characters. 


2.  Options :  As  follows: 


1)  Units : 

Enter  1  for  English  units. 
Enter  2  for  metric  units. 


2)  Normal  drag  forces: 

Enter  1  to  include  the  forces. 
Enter  0  to  neglect  the  forces. 


3)  Tangential  drag  forces: 

Enter  1  to  include  the  forces. 
Enter  0  to  neglect  the  forces. 


4)  Added  mass  forces: 

Enter  1  to  include  the  forces. 
Enter  0  to  neglect  the  forces. 
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SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


5)  Buoyancy  forces: 

Enter  1  to  include  the  forces. 

Enter  0  to  neglect  the  forces. 

6)  Gravity  forces: 

Enter  1  to  include  the  forces . 

Enter  0  to  neglect  the  forces. 

7)  Viscosities  and  fluid  mass  densities: 

Enter  0  to  use  the  default  values. 

Enter  1  to  use  substitute  values  entered  on  four 
subsequent  data  lines  as  follows: 

i)  Viscosity  of  the  first  fluid 

ii)  Mass  density  of  the  first  fluid 

iii)  Viscosity  of  the  second  fluid 

iv)  Mass  density  of  the  second  fluid 

8)  Fluid  velocities: 

Enter  0  for  zero  fluid  velocities. 

Enter  1  for  constant  fluid  velocities.  Then,  on 
the  next  two  lines,  enter  for  the  two  fluids, 
the  X  and  Z  components  of  the  velocities. 
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SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


9)  Fluid  Interface  Level: 

Enter  a  0  for  a  zero  interface  level. 
Enter  a  1  for  a  nonzero  interface  level. 
Then,  on  the  next  line,  enter  a  distance. 

3.  Number  of  links:  Enter  number. 

4.  Individual  link  data:  As  follows: 


1)  Link  number 

2)  Label:  Up  to  20  characters 

3)  Link  reference  point  label:  Up  to  20  characters 

4)  Link  number  of  adjacent  lower  numbered  link 

5)  Link  weight  density  (per  unit  length) 

6)  Link  length 

7)  Link  diameter 

8)  Initial  configuration  of  the  link:  Enter  0  the  initial 
angle  (in  degrees) 

9)  Initial  motion  of  the  link:  Enter  6  the  derivative  of  the 
initial  angle  (in  degrees/second). 

5.  Number  of  towed  spheres:  Enter  number. 
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SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


6-  Individual  towed  sphere  data:  As  follows: 


1)  Towed  sphere  number 

2)  Towed  sphere  label:  up  to  20  characters 

3)  Link  number  of  adjacent  towing  link 

4)  Towed  sphere  weight 

5)  Towed  sphere  diameter 

7 .  Number  of . towed  bodies:  Enter  number. 

8.  Individual  towed  body  data:  As  follows: 

1)  Towed  body  number 

2)  Label:  up  to  20  characters 

3)  Towed  body  reference  point  label:  up  to  20  characters 

4)  Link  number  of  adjacent  towing  link 

5)  Towed  body  weight 

6)  Towed  body  axial  length 

7)  Towed  body  volume 

8)  Towed  body  inertia 

9)  Towed  body  mass  center  position  vector  relative  to 
cable  attach  point  (X,  Z  components) 

10)  Buoyancy  center  position  vector  relative  to  the  mass  center 
of  the  towed  body  (X,  Z  components) 
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SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


11)  Initial  configuration  of  the  towed  body:  Enter  0  the 
initial  angle  (in  degrees) 

12)  Initial  motion  of  the  towed  body:  Enter  0  the  derivative 
of  0  the  initial  angle  (in  degrees/second). 


13)  Hydrodynamic  coefficients  for 

the  towed  body: 

Line  1 : 

X'  ,  X'  ,  X'  , 

and  X' 

u’ 

uu  ww  wq 

qq 

Line  2: 

X6 

and  X: 
o 

b 

s 

Line  3: 

z'«, 

Z'*,  Z'  ,  Z'  , 

and  Z' 

w 

q  uu’  uw 7 

uq 

Line  4 : 

zw|  w 

1  ’  Zqlw|  ’  Z6b’ 

and  Z 1 
o 

s 

Line  5 : 

M  , 

M'-,  M'  ,  M'  , 

and  M' 

w’ 

q  uu  uw 

uq 

Line  5 : 

M' , 

. ,  m".  , ,  m;  , 

and 

w  |w 

1  q|w|’  6b’ 

6 

s 

Motion  of  the 

mean 

ship  frame  and 

system  reference 

Enter  data  only  for  nonzero  variables  among  the  following: 


Identification  Number 
1 


Variable 

Forward  displacement  of  reference  point 
Q  relative  to  the  mean  ship  frame 

Downward  (vertical)  displacement  of 
the  reference  point  Q  relative  to  the 
mean  ship  frame 

Forward  speed  of  the  mean  ship  frame 
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SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


First  eater  the  number  of  nonzero  variables  which  are  to  be  specified. 
On  the  next  line  enter  the  variable  identification  number  in  column  1, 
and  in  column  2  enter:  1  for  a  precoded  function;  2  for  an  accelera¬ 
tion  profile. 

If  a  precoded  function  is  selected,  values  of  f^,  A,  B,  p,  and  <})  on 
the  subsequent  line.  (See  Equation  (25).) 

If  an  acceleration  profile  is  selected,  enter  the  following  data  lines: 

1)  Number  of  points  (  up  to  25) 

2)  Initial  time;  Initial  acceleration;  Initial  speed; 

Initial  displacement  (all  on  one  line). 

3)  Time;  acceleration  (one  line  for  each  remaining  paint). 

8.  Integration  Options:  Enter  on  a  single  line:  i)  starting  time; 

ii)  ending  time;  iii)  time  increment;  and  iv)  error  bound. 

9 .  Output  Options : 

1)  Enter  the  printing  time  increment 

2)  Enter  number  from  0  to  4  in  column  1  as  follows: 
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SUMMARY  OF  INPUT  DATA  REQUIREMENTS 


0  -  Print  all  data. 

1  -  Print  all  but ^tension  at  the  joints. 

2  -  Print  as  above  except  for  tension  at  the  reference  point  Q. 

3  -  Print  above  except  for  joint  kinematics. 

4  -  Print  as  above  except  for  mass  center  kinematics. 
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